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1.0  INTRODUCTION  AND  LITERATURE  SURVEY 

1.1  liitj:gdiigj:flry  Comments 

The  streamwise  evolution  of  2-D  unsteady  fluctuations  in  a 
parallel-flow  boundary  layer  is  formulated  in  this  analysis.  The 
normal  velocity  fluctuation  and  derivatives  are  specified  along  the 
y-axis,  and  the  behavior  in  the  boundary  layer  downstream  ojf  the 
y-axis  is  obtained.  Permissible  disturbances  include  the  Tollmien 
wave  and  other  eigenwaves,  but  also  include  vortical  and  irrotational 
fluctuations  in  the  freestream.  The  analysis  links  the  boundary 
conditions  with  the  amplitudes  and  phases  of  the  various  waves. 

This  boundary-value  problem _,uses  the  asymptotic  solution  of  the 
Oct-Sommerfeld  equation  by  Tsuge  (Ref,  1)  which  is  summarized  here  in 
Appendix  A.  With  ^this  asymptotic  theory,  the  Tollmien  wave  was 
calculated  by  Tsuge  and  Sakai  (Ref.  2) . 

This  report  describes  the  details  of  Ref.  3  and  is  similar  to  the 
AIAA  paper  (Ref.  4) , 

For  brevity,  the  list  of  initial-value  problems  in  time  (Section 
1.2)  and  the  list  of  boundary- value  problems  in  space  (Section  1.3) 
are  restricted  to  problems  that 

(a)  analytically  decompose  initial  or  boundary  conditions  into 
the  solutions  of  the  Orr-Sommerfeld  (or  Rayleigh)  ecpiations,  or 
those  equations  with  a  uniform  mean  flow 

(b)  include  freestream  disturbances  (i.e.  attention  is  not 
restricted  to  the  evolution  of  the  discrete  eigenmodes,  except 
for  channel  flows  where  an  infinite  set  of  discrete  eigenmodes 
is  mathematically  complete) 

The  Fourier-Laplace  solutions  of  the  Orr-Sommerfeld  equation  have  been 
documented  in  other  papers  and  reports.  A  summary  of  the  known  2-D 
and  3-D  solutions  appears  in  Refs,  5b  and  5c. 

Besides  the  studies  summarized  in  Sections  1.2  and  1.3,  other 
approaches  have  been  used  to  study  unsteady  shear  layers,  including 

(1)  Superpositions  of  eigenmodes  with  their  amplitudes  specified. 
References  using  this  approach  include  Criminale  (Ref.  6,  Chpt.  V) , 
Caster  (Ref.  7) ,  and  Mack  and  Kendall  (Ref.  8) . 

(2)  Inviscid  analyses  of  free-shear  layers  with  step-function  mean 
velocity  profiles,  e.g.  Bechert  (Refs,  9  and  10) 

(3)  Analytical  solutions  of  the  unsteady  Prandtl  boundary  layer 
equations  or  higher-order  boundary  layer  theory,  as  surveyed  by 
Teleonis  (Ref.  11) 

(4)  Lagrangian  descriptions  of  turbulence  convected  past  blunt  bodies, 
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as  described  by  the  sudden  distortion  theoriesr  e.g.  Hunt  (Ref.  12) 

(5)  Numerical  solutions  for  shear  layers  using  (a)  the  linearized, 
elliptic  momentum  equations  with  a  specified  mean  flow,  (b)  the 
parabolized  Navier-Stokes  equations,  (c)  the  Navier-Stokes  equations, 
(d)  Lagrangian  tracking  of  vortex  filaments,  and  (e)  equations 
governing  the  evolution  of  the  correlations  with  truncations  or  models 
for  the  higher  correlations.  Examples  of  these  approaches  are  cited 
in  Ref.  5c 

(6)  Compressible  boundary- value  or  initial-value  analyses  with 
freestream  disturbances  by  Tam  (Refs.  13-15) ,  Shapiro  (Ref.  31) ,  and 
Thomas  and  Lekoudis  (Ref.  17) 

(7)  Deterministic  experiments  with  vortical  freestream  disturbance  by 
Dovgal,  Kozlov,  and  Levchenko  (Ref.  18)  and  Rockwell  (Ref.  19) . 


1 .2  List  Initial-Value  Problems  In  Time  for  Patflllel  FlQWS 

(1.2a)  An  array  of  vortices  bisected  by  a  plate  for  a  general 
orientation  of  the  vortices,  including  viscosity  (Ref.  20;  Ref.  21  is 
closely  related  mathematically) 

(1.2b)  3-D  oblique  plane  waves  of  vorticity  bisected  by  a  plate, 
including  viscosity  (Ref.  22) 

(1.2c)  Free-shear  layer  with  uniform  mean  vorticity  excited  by 
vortical  fluctuations  inside  and  outside  the  shear  layer;  resonant 
and  nonresonant  excitation  (Ref,  23),  This  study  was  documented  in 
"The  temporal  response  of  a  free-shear  layer  to  vortical  disturbances" 
(1978),  an  unpublished  paper  available  from  the  author. 

(1.2d)  3-D  vorticity  in  a  boundary  layer  represented  by  a  layer  of 
uniform  mean  vorticity;  inviscid  analysis  by  P.  Durbin  "Distortion 
of  turbulence  by  a  constant-shear  layer  adjacent  to  a  wall,"  private 
communication  (1977) . 

(1.2e)  2-D  fluctuations  in  a  boundary  layer  with  two  layers  of  uniform 
mean  vorticity  in  unstable  and  stable  configurations,  inviscid 
analysis  (Refs.  24  and  25) 

(1.2f)  Boundary  layer  with  smooth  velocity  profiles,  2-D  fluctuations, 
including  viscosity  (Ref.  26) 

(1.2g)  Nonlinear  temporal  evolution  of  longitudinal  vortices  in  a 
parallel-flow  boundary  layer,  inviscid  (Ref.  20) 


1.3  List  ^  Boundary-Value  Problems  In  Space  £n£.  Parallel  aM 

Non-Parflilel 

(1.3a)  Downstream  viscous  2-D  half-plane  problem  with  arrays  of 
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vortices  (Refs.  27  and  28) 

{1,3b)  Inviscid,  uniform  mean  flow,  2-D  fluctuations  in  a 
quarter-plane  with  a  plate  (Ref.  29} 

(1.3c)  Eight  full-plane  problems  with  a  semi- inf inite  plate  (with  a 
leading  edge),  uniform  mean  flow  (Refs.  30-33).  These  solutions  are 
represented  as  a  traveling  wave  and  a  superposition  of  standing  waves 
in  Ref,  34 

{1.3d)  Full-plane  problem  with  a  finite-length  plate,  inviscid,  2-D 
fluctuations  (Ref.  30) 

{1.3e)  Inviscid  problem  with  a  mean  boundary  layer  having  two  layers 
of  constant  mean  vorticity  (Ref.  35) .  This  analysis  shows  the 
eigenmodes  excited  by  vortical  freestream  disturbances. 

(1.3f)  Viscous,  downstream  quarter-plane  problem  with  a  smooth 
velocity  profile,  2-D  fluctuations  (Refs,  3,4,36) 

(1.3g)  Viscous,  downstream  quarter-plane  problem  with  a  2-D  stationary 
wavy  wall,  2-D  fluctuations,  smooth  mean  velocity  profile  (Ref.  37) 

(1.3h}  Flow  past  a  semi- inf inite,  stationary  wavy  wall  (Ref,  38) 

(1.31)  Plow  past  a  semi-infinite  plate  with  traveling  sinusoidally 
surface  waves  that  travel  at  speeds  different  than  the  freestream 
speed,  and  galloping  surface  waves  that  travel  at  the  freestream  speed 
{Ref.  39) 

(1.3j)  Boundary-value  problem  for  a  channel  (Ref.  40} 

(1.3k}  3-D  fluctuations  in  a  viscous  boundary- value  problem  with  a 
smooth,  parallel-flow  boundary  layer  (Ref.  41) 

(1.31)  Coupling  between  an  oscillating  freestream  and  a 
Tollmien-Schlichting  wave  in  a  nonparallel  boundary  layer  (Ref.  42) 

Other  analyses  have  been  carried  out  that  link  the  initial 
conditions  with  amplitudes  of  the  instabilities.  These  studies 
include  nonresonant  and  resonant  excitation  of  buoyancy  instabilities 
(Refs.  43-45)  and  the  temporal  evolution  of  3-D  disturbances  in  an 
Ekman  boundary  layer  (Ref.  46} . 

In  the  next  section,  the  partial  differential  equation  describing 
the  evolution  of  disturbances  will  be  derived  and  integral  transformed 
to  yield  a  forced  Orr-Sommerfeld  equation.  In  Section  3,  an 
asymptotic  solution  of  this  forced  equation  is  found.  The  inverse 
Laplace  transform  is  obtained  in  Section  4.  The  formulation  is 
summarized  and  discussed  in  Section  5.  The  method  of  successive 
approximation  for  solving  the  Orr-Sommerfeld  equation  is  summarized  in 
the  appendix. 
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2.0  DERIVATION  OP  THE  FORCED  3-D  ORR-SOMMERFELD  EQUATION 
2.1  Introductory  Commentfl 

three  ^mentura  equations  for  a  constant  property  flow  are 

and  prnssur.  separated  Hto  n^n  a^l 

fluctuation  components  and  (2)  averaged  to  obtain  equations  for  tho 
^an  ^antitiea.  The  averaged  equatlona  are  subtrlSS'tJSn  Soee^f 
to  yield  equations  for  the  disturbances,  and  the 
disturbance  equations^are  linearized  for  small-amplitude  fluctuations 

and  V  =  W  =  0,  is  assumed.  DerivatJ?er  of 
equations  and  the  continuity  equation  are  combined  so 
and  two  of  the  three  velocities  arreJImiSJed?  T^e 
SrSlS^flo3“iB^^^  governing  the  evolution  of  3-D  fluctuations  in  a 


1 


dx  j  dX 


(2.1) 


where  €=  1/R  is  the  inverse  of  the  Reynolds  number. 

equation  allows  us  to  formulate  a  boundary-value  problem  in 
terms  of  solutions  of  the  homogeneous  Or r - Somme rf eld  equation.  This 
equation  results  when  you  seek  solutions  of  form 

j2.2) 

complex  function,  o<  and  y  are  the  wavenumbers  in  the  x 

“nid  titi?rJss:=;hrtrori-.‘;:^e^L«ur;eiSi?>  -  ---- 

where  D  =  d/dy  is  the  ordinary  derivative. 

primitive  variables  {such  as  velocities 
fo^^thP  t'pbs  and  pressure)  are  sometimes  more  useful  than  equations 
intuitive  variable  <^.  For  example,  equation  (2,1)  has  a 

in  tS  ?nes?r^  «  R^t.’^rindJoIt^f 

the  2-D  vorticity 

'  y  ^ 

Idf 


dt  dX  J  ? 


{2.4} 
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indicates  that  thjg  vorticity  convects  (Tj<+U?5^)r  diffuses  and 
is  produced  in  _the  shear  layer.  For  inviscid  flows  with 
irrotational  mean  flows  (f  =  0)  or  flows  with  uniform  mean  vorticity 
(fy  =  0)f  the  disturbance  vorticity  convects  with  the  mean  flow.  For 
these  flows f  the  vorticity  is  found  by  (analytically  or  numerically) 
tracking  the  vorticity  as  it  convects  downstreamf  followed  by  solving 
Poisson's  equation  for  the  disturbance  streamf unction.  Ey  this 
method,  solutions  have  been  found  which  describe  effects  of  leading 
and  trailing  edges,  lateral  edges,  trailing  vortex  sheets,  and  other 
phenomena  with  practically  iit^xjrtant  physics  and  mathematics. 
Although  links  are  possible  and  desirable  between  these  works  and  the 
Orr-Sommerfeld  equation,  direct  study  of  the  Orr-Sommerfeld  equation 
has  not  provided  the  best  insights  into  these  new  problems. 


2.2  Domains  fOI.  Boundary-Value  Problems 

Figure  2.1  shows  the  quarter-plane  region  (x  >  0,  y  >  0)  studied 
here.  Sufficient  conditions  on  the  velocity  fluctuation  are  specified 
along  the  y-axis,  and  along  the  plate  coinciding  with  the  x-axis,  so 
that  solutions  can  be  obtained  in  that  quadrant. 

The  domains  for  other  boundary-value  problems  are  also  sketched 
in  that  figure. 


2.3  Intnaral  Transforms  of  ffg  Partial  Differential  Equation 

The  coefficients  U  and  .V  in  equation  (2.1)  are  independent  of 
X,  z,  and  t.  Thus,  integral  transforms  in  those  directions  reduce  the 
partial  differential  equation  to  an  ordinary  differential  equation. 
We  focus  on  problems  with  "steady-state  oscillations"  where  our  wind 
tunnel  has  been  turned  on  for  a  long  time  and  the  conditions  along  the 
y-axis  oscillate  in  time  (but  not  necessarily  sinusoidally) .  The 
transient  response  arising  from  the  start-up  procedure  has  vanished 
from  the  region  of  interest. 


For  this  case,  equation  (2.1)  is  Fourier  transformed  in  time. 
Because  the  z-domain  extends  from  -«»<  z  <»,  a  Fourier  transform  in  z 
is  taken.  In  Chapter  3,  we  study  a  2-D  flow.  The  complex  Fourier 
transform  in  time  and  the  z-direction  is 


y,y,^)=  \  \  y( Xy  y, 


-t-  rs*  fcvt 


(2.5) 


The  signs  in  the  exponential  have  been  chosen  so  the  resultant 
transform  is  consistent  with  the  classical  Orr-Sommerfeld  equation. 
An  alternate  transform  is  the  generalized  Fourier  transform 


-  r  jre -  ^  -t-  r cut'  ~ 


(2.6) 


By  either  transform,  equation  (2.1)  reduces  to 
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upstream 

quarter-plane 


dovmstream 

quarter-plane 


half-strip  in 
y-direction 


upper  half -plane 


upstream 

half-strip 


downstream 

half-strip 


rectangle 


full-strip 
in  x-direction 


Figure  2.1  Some  domains  for  2-D  boundary-value  problems 
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oh  -  eV 

\  dX 


A  2 


dx 


(2.7) 


(2.8) 


where  the  Laplacian  reduces  to 

z. 

ax*- 

Because  we  are  interested  in  the  downstream  quadrant  in  Figure  2.1 ^  a 
one-sided  Laplace  transform  in  the  x-direction  is  taken 


r  ^  —  SX 


(2.9) 


where  s  is  the  complex  Laplace  parameter.  With  this  transformation, 
the  following  equation  results 

.2  ^ 

^  (2.10) 


When  s  =  i><is  introduced,  where  o(  is  the  complex  x-wavenumber ,  then 
the  bracketed  term  on  the  left-hand  side  is  the  classical 
Orr-Sommerf eld  operator  in  equation  (2.3) .  The  right-hand  side  is  the 
forcing  function 

(2.11) 

where  v*^s  the  normal  velocity  (or  its  derivatives)  along  the  y-axis. 
For  a  2-D  flow  with  IT  =  0 ,  the  forcing  function  is 

(2.12) 

To  introduce  definiteness  in  our  analyses,  the  forcing  function  is 
Fourier  transformed  and  represented  as 


(2.13) 


Thus,  the  forced  Orr-Sommerf eld  equation  for  a  2-D  flow  studied  in  the 
next  chapter  is 

2, 


S  (2.141 

where  the  Laplace  parameter  and  the  x-wavenumber  are  related  by  s=ip/. 
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3.0  SOLUTION  OP  THE  TWO-DIMENSIONAL  FORCED  ORR-SOHMERFELD  EQUATION 

3.1  Particulac  Integral  ani  Asymptotic  Form 

The  objective  now  is  to  obtain  the  particular  solution  of  the 
forced  Orr-Sommerfeld  equation 

where  the  forcing  function  is  prescribed  at  x  =  0  as  the  Fourier 
component 

The  particular  solution  is  obtained  by  the  method  of  undetermined 
coefficients.  For  the  Orr-Sonmierfeld  equation,  the  solution  was  given 
by  Gustavsson  (Ref.  26)  as 


(3.3) 


where  W  is  the  Wronskian 


tV  = 


< 


(3.4) 


and  4’  3  =  1,2, 3, 4  is  the  cofactor  of with  the  sign  reversed. 

for  j  =  1,2, 3, 4  are  four  independent  solutions  of  the  Orr-Sommerfeld 
equation. 


Because  the  Orr-Sommerfeld  equation  lacks  the  next-to-highest 
derivative  (the  third  derivative),  then  the  Wronskian  is  constant. 
This  meeuis  that  the  Wronskian  only  needs  to  be  evaluated  at  one  point, 

y  =  y, . 

The  wall  conditions  and  the  asymptotic  solutions  for  large  y  are 
listed  in  Table  1.  The  four  solutions  are  derived  in  Appendix  A. 
These  asymptotic  solutions  are  based  on  a  new  scaling,  and  are  for 
flows  where  0<y<®or-a><y  <»  . 

Substituting  the  asymptotic  expressions  for  the  tpB,  then  the 
Wronskian  is 

,,  -2  ,  i/z 

(4^=  4^'*''^  v^-^j  where  ~  (3.5) 


12 


AEDCTR83-9 


Characteristics 

Table  1. 

of  the  four  independent, 

asynptotic 

solutions 

of  the  Orr-Somirerfeld  equation 

Solution 

Value  at 

y=0  Asymptotic  form  for  y»y 

0 

— /  ^ 

(1-c)  sinho<(y-y) 

0 

(l-c}cosho((y“y) 

1 

^3  exp[-A(y-y)  ] 

0 

exp[H(yfi'{y-y)  ] 

where  = 

[«v®+iR( 

io)  ]  and  ^  =  ^3(9) 

with  the  real  part  of/**  chosen  as  positive.  Alsor  the  wall  value  of 
the  ^rticular  solution  (equation  3,3}  must  satisfy  the  impermeability 
condition,  <^(0)  =  0,  by  choosing  the  lower  bound  of  the  integral  as 
follows 

y  ®  y 

^  ^  i  a,  (3.6) 


The  asymptotic  particular  solution  is  obtained  with  the  solutions 
in  Table  1.  The  result  is  *  -* 

(»<(o^+<-j8)  t-  ' 


(3.7) 


This  solution  is  used  later  when  an  outer  boundary  condition  is 
satisfied.  Note  that  equation  (3.7)  includes  terms  which  grow 
exponentially  like  exp(«/y)  .  These  terras  are  cancelled  later  so  that 
the  complete  solution  is  well  behaved.  Equation  (3.7)  does  not 
include  terms  like  exp(i^y). 


3.2  COfliBlementary  Solution  and  tllfi  Assembly  JJifi  General  Solution 

The  general  solution  of  the  forced  Orr-Sommerfeld  equation  (3.1) 
is 
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J^f 

subject  to  the  boundary  conditions 


=  0  {inpermeability) 
=  0  (no-slip) 
and  if  is  bounded  as  y 


at  y  =  0 


(3.8) 

(3.9) 

(3.10) 

(3.11) 


Because  four  of  the  five  solutions  satisfy  in^aermeability 

then  for  the  general  solution  (3.8)  to  satisfy  impermeability,  then 
the  coefficient  of  the  fifth  term  which  does  not  satisfy 
impermeability  must  vanish 

<^3  =  (3.12) 

Also,  because  only  ^  grows  like  exp^y)  as  y  then  the 

boundedness  condition  (3,11)  requires  that  coefficient  to  vanish  also 


and 


Sl.  +  c  (f-c)-i- ~ —  r 


'^9 


(3.13) 

determined  by  the 
n  3.11)  conditions 

no-slip 

O 

II 

(3.14) 

-V-  - ; —  /  =  o 

*S<+  (  ^  *' 

(3.15) 

the  coefficient  of  exp(a(y) 
Note  that^does  not  vanish  as 

y  ^at,  but  oscillates  as  exp(±i|Sy). 

These  two  equations  for  the  two  constants  are  solved  to  give 

^)] 

l+eB  yJ 


o/-  e 


-'F 


af-t-njS 


(3.16) 
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o  ^ 

where  ^'(o)  =  ^/(o)  (  4,  J  '^'&  ^  \  (3.17) 

.7  00 


if 


A 

y 


<f^/C^)  %(o) 

{ck-coY 


(3.18) 


Note  that  4  =  0  is  identical  with  the  eigenvalue  condition  for  the 
discrete  modes,  as  seen  from 


<p'^  C,  <f>/  <^^=.  c? 


y 


and 


^  —  e  ^  C,  - €■  =  C7 


y 


ao 


'  (3.19) 

The  impermeability  condition  is  automatically  satisfied  by  and  <fi^  . 


A  more  precise  look  at  the  behavior  of  reveals  that  must  be 
rescaled  in  the  non-asyitptotic  region.  The  behavior  in  the 
non-asyitptotic  region  and  the  asynptotic  (outer)  region  are 

/-  c 


where  Bois  a  constant  depending  only  on  c  =<J6'/4i;  it  is  not  affected 
by  rescaling.  Thus  must  be  rescaled  so  that 


(3.20) 


for  the  non-asynptotic  expression  for  <5^  . 

The  determinant  A  of  equation  (3,19)  has  multiple  zeros 
corresponding  to  the  eigenvalues  of  the  discrete  modes.  The 
Tollmien-Schlichting  mode  is  the  lowest  or  fundamental  mode.  Then 
1/dis  expanded  as  the  partial  fractions 


J_ 

A 


E 


A/ 


fN) 


with 


/ 


(3.21) 
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Using  i^as  rescaled  by  equation  (3.20),  then 


^  A  —  ^ 


oi 


<P,  Co) 

^oi 

«✓ 


^^Co) 

1  + 


/  zC^-ee}) 


CN) 


Because  4  =  0  at  ,  then  h^is 


<P/0)/^  ot  ^  <P'(0)A^ 

o(  (oi-acf)^ 


+  o( 


4^/(0)  %^(o) 

/  z(=i~ca') 


-\ 


(3.22) 


3 ,3  Suimiiarv  ^  the  solutions  in  (ytawst^  Space 

The  final  expression  for  the  general  solution  of  the 
Orr-Sonrmerfeld  equation  is  ^ 

I 

(3.23) 

The  new_<^as  rescaled  in  equation  (3.20)  is  used.  Because  $  is  linear 
in  ^  ,  it  is  rescaled  by  the  same  factor 


'  ana  ^  (3.24) 

Far  from  the  wall  (y>>y) ,  an  asymptotic  general  solution  results  when 
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the  asymptotic  expression  (3.7)  is  substituted  into  equation  (3.23) . 
The  result  is 


~^Cy'y) 


(E^-)  l 


I—  %'(o) 


.  X  e-  ^  ro;  ]  \r.. 


<?o<  —  £<J  ^  ^  f  js  •<1 

o(  -*  f  S  Z  /  -3  oc\ 
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4,0  INVERSE  LAPLACE  TRANSFORM 

4.1  General  Expression  tha  Contributions  from  Poles  Branch 

Lines 


In  Chapter  3,  the  general  solution  of  the  forced  Orr-Sommerfeld 
equation  in  transform  space  was  obtained.  To  express  the  fluctuation 
velocity  as  a  function  of  x  and  y,  the  inverse  Laplace  transform  is 
necessary 

-  /  ^ 

K  =»  —  \  </s 

r  (4.1) 

where/* consists  of  the  contours  s  =  s<,+  iA  and  a  semi-circle  of 
infinite  radius  in  the  domain  of  Real(s)  <  0,  The  real  value  s  is 
chosen  so  that  the  poles  of^are  located  on  the  left  side  of  the  line 
So+  iA  rWith  -a><A<<»,  Figure  4.1  shows  this  contour,  the  poles,  and 
the  branch  lines. 


The  contribution  from  the  jth-pole,  s  =  s : ,  is  calculated  from 
the  Cauchy  theorem 


r 


CJ) 


(4.2) 


If  <*^is  not  single-valued  in  the  s-plane,  as  caused  by  fractional 
powers  of  s,  then  a  Riemann  cut  is  required.  Because  the  integration 
contour/* cannot  pass  through  the  Riemann  cut,  it  must  detour  along  the 
branch  cuts,  Y  ,  or  along  another  path  which  does  not  cross  those 
cuts.  This  rule  gives  us  the  following  expression 


O’j 


(4.3) 


The  contributions  from  the  poles  and  branch  line  is  obtained  in  the 
following  sections. 


4.2  Interior  Solution 

4.2.1  Poles  Iiom  i±&  Discrete  Eigenmodes 

Bqn.  (3.23)  for  the  interior  solution  shows  that  there  are^  N+2 

poles  corresponding  to  s  =io/j^*'and  s  =:t^if  the  forcing  function  F’is  a 

pole-free  analytic  function.  The  equation  also  shows,  for to  have  a 

positive  real  part,  that  branch  lines  are  required  because  of  the 

fractional  power  in  eqn.  (3.5).  Because  F*  is  open  to  free  choice, 
poles  or  branch  lines  can  be  generated  or  eliminated  by  it. 

The  poles  s  =  i°/,  of  eqn,  (3.23)  are  for  the  discrete  eigenmodes. 
According  to  the  Cauchy  theorem  (4.2),  we  have 
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o^  $'Co) 
JZ/s//r 


•^+r/ 


r  jCo) 


+  <’■><' 


F  e 


f 


(4.4) 

Each  mode  satisfies  the  impermeability  and  no-slip  conditions 
independent  of  each  other  mode  and  independent  of  other  contributions. 
These  conditions  are 


^  ^  and  ^ 

The  first  randition  is  easily  checked,  <^(0)  =<?{0)  =  0.  The 
condition  is  fulfilled  because  of  the  eigenvalue  condition 


^/c°) 


which  vanishes  for  «*/  =o< 


<n) 


(4.5) 

second 


(4.6) 


4.2.2  Poles  icam  tiie  standing  Wave  Modes 

The  contribution  from  the  pole  s  =  -|3=  Sj(or  =  ifi)  is  also 
obtained  in  a  straightforward  manner.  The  result  is  ' 


(4.7) 

This  mode  is  a  standing  wave  which  decays  as  exp(-8x).  It  satisfies 
impermeability  and  no-slip.  This  mode  was  studied  by  Rogler  and 
Reshotko  in  the  freestream  (Ref.  27)  and  appeared  in  another 
boundary-value  problem  in  a  quarter-plane  (Ref.  29) .  Solutions  of  the 
Orr-Sommerfeld  equation  for  boundary  layers  were  calculated  by  Rogler 
and  Tsuge  (Ref.  5a)  and  Rogler  (Refs,  5b— c)  for  this  decaying  wave  and 
the  corresponding  growing  wave.  This  growing  wave  is  discussed  next. 

The  contribution  from  the  exponentially  growing  standing  wave. 


f 
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s  (or  **<  ~  -is)  can  be  calculated  in  the  same  manner.  However,  this 
mode  is  suppress^  in  this  quarter-plane  problem  where  x  .  To 
exclude  this  mode,  the  forcing  function  must  have  the  form 

^  (4.8) 

where  m  is  an  integer.  Physically  realizable  situations  exist  where 
this  growing  mode  is  acceptable,  such  as  the  upstream  influence  of  a 
trailing  edge  or  another  downstream  boundary  condition.  The 
possibility  of  modeling  such  effects  in  a  a  problem  with  0<x<*is  not 
considered  here.  This  growing  mode  can  appear  in  boundary-value 
problems  where  the  x-domain  is  finite  with  0<x<L,  or  semi-infinite 
with  -®<x<0,  as  Figure  2.1  shows. 


4.2.3  Contributions  Lins.  yielding  e  Continuous 

Saectrum 


As  has  been  defined  in  eqn.  (3.5),  the  function 


^  S^-f-  /^s  —  f  ) 


(4.9) 


is  not  single-valued  in  the  s-plane.  A  branch  cut  is  required  which 
makes  ^  unique  and  maps  the  whole  s-plane  onto 


(4.10) 


The  contribution  from  this  branch  line  is  an  integral  over  the 
continuous  (vortical)  mode.  Condition  (4.10)  requires  the  branch  cut 
to  be  a  special  curve.  To  simplify  the  task  ahead,  we  let 


where 


[(S-  SgXSj- 

Z  \  '/^ 

s"=  ('T  - 

®  2 


(4.11) 

(4.12a) 


(4.12b) 


and  define 
Then 


(T  -  S  -  S3 

[^^(53-  53-0*;]^ 


(4.13) 

(4.14) 


Now  consider  a  curve  along  which  Real^)=0  (4.15) 

Along  that  curve  a/'^  (T'  ayy(s'-  s^-cr')—  ±  (4.16) 

The  two  vectors  (T  and  -  s  ~  <T  and  their  arguments  ^and^'^are  shown 
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in  Figure  4.2.  The  locus  of  the  curve  is  obtained  from  the  following 
relationship  by  taking  the  tangent  of  (4.16) 


£  ^  (Sj"-  s,).  -  ^  ^ 

(  Sj'  - 


(4.17) 


This  branch  line  forms  part  of  a  hyperbola,  labeled  in  Figure  4.1, 
along  where  is  purely  imaginary. 

r  €r  ^ 


+  eT- 


t/2 


s/]  (s^  +  7-; 


/^  = 


Sf.+zr 


T  »• 


Sr  = 


S  •  = 


(4.1B) 


Along  its  upper  side,  rmag(/^)  >  0  and  along  its  lower  side, 
lmaq[/f)  <  0.  The  whole  s-plane  is  mapped  onto  a  half  plane 
Real{^)  >  0  by  this  branch  cut.  The  contribution  from  the  branch  cut 
is 


,  -C°\  c/^Cs^-^zr-)  j[s3-r(l-  | 

(4.19) 


evaluated  at  «V  =  _fS  ] 

^  '  5^*T 


and  where  g_is 


5,0;= 


,2  ^'(O) 


(o<-6oy 


.J _ 


(4.20) 
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Contribution r  13/  satisfies  the  imperroeability  and  no-slip  conditions 

9^(0)  ^  CoJ  =  (o) ^(oj  ^  oj 

V  / '(c)  1  .  (a)  )  ^'Co) 


$'(o) 

§^(o) 


o 


z  -  2 

For  large  Reynolds  numbers,  the  terms  (s^-/s^)  0(R  )  are  ignored. 

The  expression  simplifies  to 


=^- 


S,)( 


zr 


\ 


t'/®  f  r+  } 


T/z- 


‘^=~/s^*rr^ 


iJc 
(4.21) 


rJ 


The  main  part  of  the  x-dependence  of  I3  comes  from  the  term  exp(Sax). 
Thus,  the  contribution  from  the  branch  line  is  approximately  ^ 


4 


(4.22) 


where  Sjis  given  in  eqn,  (4.12)  .  This  approximate  solution  propagates 
with  the  freestream  velocity  and  decays  slowly  in  the  x-direction. 
For  X  »  1,  the  x-dependence  of  the  integral  shows  a  universal 
feature.  The  contribution  from  the  integrand  is  limited  to  a  region 
close  to  T  “  0.  Transforming  the  variable  from  r  to  t  such  that 

■^r 

we  have 


SgX 


ao  ^  ^ 


which  shows  the  enhancement  of  the  spatial  decay  by  the  factor  x 


(4.23) 

-'/2 
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4.3  Exterior  (Asymptotic)  Solution 

4.3.1  Poles  from  thg  Dl^rete  Eigenmodes 

The  procedure  used  in  Section  4.2,1  also  applies  to  the 
asymptotic  solution f valid  for  y  »  y  as  given  in  eqn,  (3.25) 


^  e 


I  Cw- ^ 


■IS  J  (<W-  UJ)  ^  ^  ( 

(4.24) 


imJ 


tN) 


This  expression  shows  that  this  mode  decays  as  exp (-■>/,  y)  toward  the 
freestream  as  it  should. 


4.3.2  thSi  Standing  Waves 


The  asymptotic  solution  for  the  pole  at  s  =  -p 


A 

/S'. 


'  .jc  '  •  -  •  \_z_— 

—  [4^  ^ 1 -*•  - , 


(4.25) 


Because  I^has  the  form  exp(~px+ipy) ,  far  from  the  boundary  layer  edge^ 
this  mode  is  oscillatory  in  y  and  is  irrotational .  This  agrees  with 
Refs.  5a-c. 


A  similar  expression  can  be  obtained  for  the  pole  at  s  =  6 , 
corresponding  to  the  e  r^mode.  If  condition  (4.8)  applies,  then  this 
mode  is  suppressed  in  the  asymptotic  and  general  solutions. 


4.3.3  lianch  Lin£  Yielding  ^  Continuous  Spectrum 

The  contribution  from  the  branch  line  for  the  asymptotic  solution 
is  obtained  from  eqn,  (4.19)  with  the  function  gjreplaced  by 


ra  \  5  \_o<A^  ) 


ifsy/n 

evaluated  at 


o<'=  -  /5,+  cr(\-,ll£. - ) 

^  ^  +  r  y 


(4.26) 


Because  s_  «  haJ-^o  /R  from  eqn.  (4,12),  this  mode  varies  with  y  as 
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\  '^  y  (4.27) 

where  T  is  the  variable  introduced  in  eqn.  (4,18). 


4.3.4  h  £Qle  Rfipxfisenting  a.  Solution  Periodic  in  siowiy  Pf^pavinfr 

The  asymptotic  expression  (3.25)  reveals  a  pole  at 

which  has  no  counterpart  in  the  interior  solution.  It  is  located  on 
the  hyperbola  branch  line,  and  is  the  mode  forced  by  the  periodicity 
at  the  initial  line,  x  =  0,  and  continues  to  y  ■‘co  with  wavenumber  p . 

Eqn.  (4.28)  is  rewritten  as 

=  -  =  O 

The  contribution  from  the  pole  s  =  s^is 


(4.29) 


- 


^  -at 

From  eqn.  (4.29) ,  the  asyn^Jtotic  behavior  is 


^  4 


ex>o 


[  ^ed  7*)±  - ^  r  X  j 


(4.30) 


(4.31) 


This  solution  is  periodic  in  y,  travels  downstream  with  the  freestream 
velocity,  and  slowly  decays  exponentially. 


The  same  procedure  for  the  pole  at  s  =  s  leads  to  a  solution 
which  grows  like  exp(s^x).  Although  this  mode  can  be  physically 
realized  as  the  upstream  influence  of  a  disturbance,  it  is  excluded  in 
this  quarter-plane  problem  by  factorizing  the  forcing  function  as 


/? 


(4.32) 


where  n  is  an  integer. 
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4.3.5  The  cases  «/=  H  and  ^ 

The  exterior  solution  (3.25)  has  factors  ^  and  — ^in  the 
denominator,  and  thus  have  pole-like  characteristics.  The  first  one, 
however,  is  ruled  out  because  it  gives  a  solution  which  approaches  a 
nonvanishing  constant  a  .  ^ 


rw,. 


li  O 


''y^co 

This  reflects  the  fact  that  F  should  have  a  factor  of 

cV(S 


(4.33) 


in  realistic  situations,  as  inferred  from  the  form  of  eqn.  (3,1) . 

The  other  case  of  -cO,  or  c  =  1,  should  be  treated  separately. 
It  is  not  coveted  by  the  present  analysis.  In  that  case,  two 
solutions  c^and  ^  need  to  be  replaced  by  yexp(+-<y)  which  leads  to  a 
different  expression  for  the  Wronskian. 
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5.0  SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 


5.1  Summary 


■nie  boundary-value  problem  in  space  for  the  evolution  of 
two-dimensional  fluctuations  in  a  boundary  layer  is  formulated.  With 
velocity  fluctuations  specified  along  the  y-axis,  and  impermeability 
and  no-slip  imposed  along  the  plate  (which  lies  along  the  x-axis) ,  we 
seek  the  small-amplitude  solutions  in  the  quarter-plane  downstream  of 
the  y-axis.  These  velocity  fluctuations  are  functions  of  "exponential 
order"  that  can  be  Laplace  transformed. 


The  differential  equation 


O 


is  Fourier  transformed  in  time  and  Laplace  transformed 
x-direction.  The  resultant  forced  Orr-Sommerfeld  equation  is 


in  the 


(5.1a) 


where  the  forcing  function  is 


/I  (o) 

y  ■+ 


\  s 


^  {5.1b) 


The  forcing  function  is  Fourier  transformed  in  the  y-direction  and 
represented  as 


(5.2) 


The  particular  integral  of  the  forced  equation  is  obtained  by  the 
method  of  variation  of  parameters. 

The  four  independent  solutions  of  the  homogeneous  equation  are 
obtained  as  asymptotic  solutions  of  the  Orr-Sommerfeld  equation  with 
an  asymptotic  expansion  ino('/R  whereof  is  the  x-wavenumber  and  R  is  the 
Reynolds  number.  The  smallness  of  this  parameter  under  ordinary 
conditions  assures  rapid  convergence.  These  four  independent 
solutions  have  unique  properties  exploited  in  the  present  solution. 

The  complementary  and  particular  integrals  are  combined  into  a 
general  solution  which  satisfies  impermeability  and  no-slip  at  the 
wall  and  boundedness  far-away. 

The  inverse  Laplace  transform  is  found  via,  the  Cauchy  theorem 
with  a  contour  integration  taken  to  the  right  of  all  poles,  along  a 

semi-circle  of  infinite  radius,  and  with  indentations  about  the  branch 
lines. 
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The  left  and  right  handed  branches  of  one  Riemann  cut  are 
hyperbolas.  Integration  about  the  left  handed  cut  yields  the 
contribution  from  the  slowly-decaying  vortical  fluctuations  in  the 
f reestream. 

Since  contributions  from  the  right-handed  leg  are  physically 
unreasonable  with  a  plate  which  extends  to  infinity  downstream, 
conditions  are  imposed  on  the  forcing  function  to  eliminate  those 
contributions. 

Contributions  from  the  isolated  poles  for  the  discrete  modes  are 
obtained  from  the  Cauchy  theorem.  The  amplitude  of  the  Tollmien  wave 
is  found  as  the  first  in  this  series. 

The  amplitude  of  the  decaying  standing  wave  with  another  pole  is 
similarly  found. 

The  exponentially-growing  standing  wave  is  excluded  by 
introducing  a  condition  on  the  forcing  function.  This  contribution  is 
excluded  for  the  same  reason  as  given  above  for  excluding 
contributions  from  the  right-handed  branch  line. 

The  branch  lines  for  the  interior  solution  include  the  left  and 
right-handed  hyperbola  branch  cuts.  The  fluctuation  associated  with 
the  right-handed  branch  line  was  calculated  by  Rogler  (Ref,  47)  as  a 
solution  of  the  Orr-Sommerfeld  equation. 

The  poles  in  the  interior  solution  include  the  N  discrete  modes, 
the  two  modes  corresponding  to  the  exp(ifix)  standing  wave  modes,  and 
poles  along  the  hyperbola  branch  lines,  ' 

The  poles  in  the  asymptotic  outer  solution  include  the  tails  of 
the  N  discrete  modes,  the  oscillatory  solutions  corresponding  to 
,  and  the  oscillatory  solutions  corresponding  to  the  standing 
'waves  with  s  =  -HicV. 


5.2  Procedure  £&r.  Calculating  Jtiie  initial  Amplitude  Tollmien 

Wave 

(1)  The  procedure  for  calculating  the  four  independent  solutions  of 

the  homogeneous  Orr-Sommerfeld  equation,  ^  f  and  is 

outlined  in  Appendix  A  and  Ref.  1.  These  <^s  are  solutions  of  the 
Orr-Sommerfeld  equation,  as  asymptotic  solutions  or  as  numerical 
solutions  with  edge  conditions  defined  by  the  solutions  of  Table  1. 
The  derivatives  and  Dj^are  evaluated  at  y=0,  and  the  ^s  up  to  the 
2nd  derivative  are  stored  for  0<y<y. 

(2)  The  derivatives  a  (di^/dy}/d^ and  a {d<;^ /dy) /So^are  evaluated  at  y  =  0 
and  for  ssV,  corresponding  to  the  Tollmien  wave  and  stored. 

(3)  h^is  evaluated  from  equation  (3.23) . 
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(4)  |(0)  is  calculated  from  equation  (3.17) . 

(5)  l/^\o,y)  is  evaluated  from  equation  (4.4)  for  the  amplitude  of  the 
Tollmien  wave  at  x=0. 


5.3  Conclusions  Recommendations 

A  boundary-value  problem  is  formulated  which  illustrates  how 
initial  fluctuations  are  represented  as  a  superposition  of  its 
(spatial)  Fourier-Laplace  solutions.  Figure  4.1  serves  as  a  roadmap 
for  the  five  families  of  fluctuations  which  can  appear  in  a  boundary 
layer.  The  solution  is  complete  in  the  sense  that  a  superposition  of 
these  fluctuations  completely  describes  the  fDowfield  for  an 
incompressible,  viscous,  small-amplitude,  2-D,  unsteady  flow  in  a 
parallel-flow  boundary  layer  over  a  flat  plate. 

Since  this  formulation  in  1980,  these  five  classes  were 
calculated  as  solutions  of  the  Orr-Sommerfeld  equation.  Their 
streamlines  are  plotted  in  Figures  B.l-B,5  in  Appendix  B  of  this 
report.  A  description  of  each  fluctuation  appears  beside  each  figure. 
These  descriptions  provide  an  overview  of  the  five  families  that  serve 
as  basic  building  blocks  for  fluctuations  in  boundary  layers.  The 
five  are; 

1.  The  diS£r£t£  eioenmodes.  where  the  Tollmien  (or 
Tollmien-Schlichting)  wave  is  the  fundamental  mode  of  boundary 
layer  instability.  These  fluctuations  vanish  far  from  the 
boundary  layer . 

2.  The  slowly-decaying  vortical  fluctuations,  with  solutions 
described  and  calculated  by  Rogler  and  Reshotko  (Ref.  48)  and 
Salwen  and  Grosch  (Refs.  49,50).  Far  from  the  boundary  layer, 
these  fluctuations  are  rotational. 

3.  The  explosively-growing  vortical  fluctuations  representing 
upstream  diffusion  of  vorticity  (Refs.  27,47).  Par  from  the 
boundary  layer,  these  fluctuations  are  rotational. 

4.  The  exponentifllly^dficaying  standing  described  in 
Refs,  5a,b,c,  In  the  freestream,  these  waves  are  irrotational. 
With  the  help  of  Mr.  Arnie  Rosner,  a  16  mm  animated  movie  was 
produced  which  illustrates  these  waves.  The  movie  shows  that 
these  waves  travel  in  the  direction  perpendicular  to  the  wall. 
Ref,  34  shows  that  a  spectrum  of  these  standing  waves  appear 
downstream  of  the  leading  edge  when  freestream  disturbances 
interact  with  that  leading  edge.  The  waves  can  be  generated  by 
surface  waviness  which  is  steady  (Ref.  38)  or  unsteady  (Ref.  39) 
if  a  leading  edge  exists. 

5.  The  exponent ially-growing  standing  waves  also  described  in 
Refs.  5a,b,c,  Although  these  waves  are  physically  realizable  in 
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some  problems^  they  ace  excluded  in  this  problem  with  a 
plate  that  extends  to  infinity  downstream*  Far  from  the  boundary 
layer,  these  growing  standing  waves  are  irrotational. 

In  future  studies,  the  theory  could  be  applied  to  identify  wtet 
properties  of  the  initial  disturbances  are  responsible  for  exciting 
the  Tollmien  waves,  and  the  role  of  the  mean  velocity  profile  {a.B 
influenced  by  the  mean  pressure  gradient  and  wall  roughness} *  Further 
study  is  also  required  to  account  for  the  leading  and  trailing  edges, 
3-D  disturbances,  and  the  growth  of  the  boundary  layer. 

The  3— D  counterpart  of  this  analysis  was  described  in  Ref.  41, 
but  additional  cuialysis  and  calculations  are  required  to  complete  that 
study. 
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Appendix  A 

METHOD  OF  SUCCESSIVE  APPROXIMATION 
FOR  SOLVING  THE  ORR-SOHMERFELD  EQUATION  USING  NEW  SCALING 


A.l  Iiit,rgdaci:.prY  .Comments 


This  appendix  summarizes  a  method  used 
Orr-Sommerfeld  equation 

to 

solve 

the 

(  / 

0 

1 

(A.l) 

where  D  =  d/dy. 

R  =  the  Reynolds  number. 

oi  and 

are 

the 

dimensionless  wavenumbers  in  the  streamwise  and  spanwise  directions, 
respectively,  and  A^  + 

Heisenberg  (Ref.  52)  observed  that  the  Orr-Sommerfeld  equation 
can  be  solved  analytically  in  the  inviscid  limit 


o*'  ^ 


with  a  small  wavenumber 


(A.2a) 

(A.2b) 


Tsuge  (Ref.  1)  recognized  that  the  reduced  (but  still  viscous) 
Orr-Sommerfeld  equation 


1  £  -  ^  ^ 


can  be  integrated  once  to  yield  a  third-order  equation 

.3 


(A.3) 


(A.4) 


This  equation  is  the  zeroth-order  equation  associated  with  the 
following  expansion 

.a  .  .  /  .2  ryX 


<f(9)  =  ^  (^)  (7^ 


-  r ' 
> 


with  y  =  and  ^ 

The  functions  (N=0,l,2,. ..)  obey  the  following  set  of  equations 


(A.5) 

(A.6a,b) 


=  O 


.  z  "'V 


fo)  (O 

where  the  operators  L  and  L  are  defined  by 


(A.75 
(A.8) 
(n>2)  (A. 9) 

(A.IO) 
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2  -- y  cj) 

Equation  (A«4)  is  the  zeroth-order  equation  in  this  set. 


(A. 11) 


A.2  Efiiu:  independent  solutions  zeroth-order  equation 


Equation  {A, 3  or  A.8) ,  or  its  integrated  version  (A. 4) 

Y  -  e  ( 0~  C  )Y  e  (yy  ^  c'C^  (A. 12) 


has  four  independent  solutions,^  ,/=l-4.  Solutions  and  (J^both  tend 
to  slowly  varying  functions,  whereas  ^ and  exhibit  rapid  decay  or 
growth.  The  solutions  are  found  in  the  order  in  a  form 
where  each  solution  requires  at  most  knowledge  of  the  from  earlier 
solutions.  ^  and  ^are  found  exactly  in  terms  of  ^  and  . 


Sglutifln 


The  homogeneous  solution  of  equation  (A.  4)  is 
Because  of  its  prospective  exponential  decay,  the 
assumed  of  form  y 

-  ex/’  \>^  c/y 
o  ' 


found  first, 
solution  is 


(A. 13) 


With  this  transformation  of  variables, 
equation  (A,12)  is 

A  “*■  3XX  +■  A  X 


the  homogeneous  form  of 


(  (J  ■=  o 


(A.14) 


or  alternately  and  more  conveniently 

A  =  /^  -  ^  T  cj 

^  ~  ~5Xy^->-^e'Cc/-c)X-~  ^  (A.15a,b) 


Since  we  are  seeking  a  solution  of,  .(A.15a,b)  which  d*»nayR 
exponentially  in  the  asymptotic  limit  (  A  0) ,  the  root  for  A  with 
negative  real  part  is  ,  y  , 

A  ~  -  e  {a-c) 


/<  ~  c  ( t/-  C ) 


{A,16a,b) 


In  solving  equation  (A.15a,b>  with  asymptotic  conditions  (A.16a,b) ,  it 
is  useful  to  introduce  a  variable  transformation 


3(  (a-c)~  (A.17a,b) 
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and  work  with  the  equations  in  the  new  dependent  variables  (V,W) 

l/-^A  U'-Z/  ^  ^  O 

bV  ZXW-Sc'  ^  O  (A.18a,b) 


subject  to  the  asymptotic  conditions 


O  os  ^  oo 

’  A, 

The  integration  may  be  started  at  a  point  ^  {>>1) 

values  .  , 


\/  -S’  tZc  A 

3A  ^^ZA 


(A.19a,b) 
with  the  initial 

(A.20a,b) 


and  with  /\  given  by  equation  (A, 16a)  are  moderately  small.  The 
quadrature  marches  inward  from  this  point,  where  the  nonlinearity  of 
equation  (A. 18)  appears  in  the  form 


C  Zt  ^ e’ 


f/'Z 


(A. 21) 


A 

and  is  taken  into  account  for  ^  <  p'. 

This  method  of  numerical  integration  of  the  original  equation 
(A. 15)  works  successfully  only  for  the  decaying  solution  A 

formally  identical  procedure  for  the  growing  solution  ^  suffers  from 
numerical  instabilities.  The  two  solutions  are  different  because  the 
pointy  =07  is  a  saddle  singularity  of  equation  {A. 15)  for  the  case 
treats,  but  it  is  a  nodal  singularity  from  which  infinitely  many 
solutions  emerge  for  ^ . 


A  straightforward  calculation  from  (A, 14,  A. 21,  and  A. 20)  leads 
to  the  solution 


The  classical  counterpart  of  given  in  Ref.  52  is 
T-  -r 


OD  «> 

where  8,^3  is  the  Hankel  function  of  the  first  kind  and 
defined  by 


(A. 23) 
T  and  are 


7'  - 


,1/3 


{A. 24) 
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with  subscript  c  signifying  the  value  at  the  critical  layer.  For 
T»\f  (A. 23}  takes  the  form 


If  we  employ  the  linear  approximation  (A, 24)  for  U-c  at  the  critical 
point,  solution  (A. 22)  reduces  to  (A. 25) . 


Solution 

After  one  of  the  solutions  of  the  third  order  equation  (A, 13)  is 
found,  the  difficulties  in  solving  for  the  remaining  solutions  are 
considerably  reduced  because  of  the  following  theorem;  If  n 
independent  solutions  of  the  mth  order  linear  differential  equation 
are  at  heuid,  then  the  (n+l)th  solution  is  obtained  by  solving  an 
equation  of  (m-n)th  order.  To  use  this  theorem  for  the  second 
solution  ^  requires  that  an  equation  of  the  order  3  -  1  =  2  be  derived 
by  putting 


and  by  substituting  into  equation  CA.13) ,  This,  in 
equivalent  to  claiming  that  if  G  is  assumed  in  the  form 

K 

A  +  Sj  ^ 


7 


(A.27) 
turn,  is 

(A. 28) 


then  the  transformed  variable  S  obeys  a  first  order  nonlinear 
equation.  A  simple  calculation  actually  confirms  the  assertion, 
yielding  a  Riccati  equation 


5*^  +  AS  -  W 


(A.29) 


where  A  and  W  have  been  solved  from  eqns.  (A. 8)  ^nd  (A.ll) , 
respectively.  Of  the  two  asyn^jtotic  roots  of  U.29)  for  S-O,  the  one 
that  vemishes  as  ^  -♦ » ,  namely 

5  =.  w/>y  (A. 30) 

is  the  correct  choice,  since  only  this  root  insures  that  ^  meets  the 
need  for  having  a  slow  variation  as  ^  In  fact,  then 

t  (Z 


(7  ' 


with 


A 


(£y-c) 


-3/z 


ex/>  - 


? 

^  S 


rr 


»/ 


(A. 31) 


r 
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Because  the  integrand  in  (A. 31)  is  a  rapidly  growing  function, 
the  portion  near  the  upper  bound  contributes  to  the  integral, 
the  integral  is  replaced  with 


only 

Then 


It  is  interesting  to  compare  of 
counterpart 


r«s;. 


/V- 


(A. 31) 
-./5 


with 


its  classical 


^ ^  4^)  (A.32) 

which  represent  the  inner  (viscous)  and  outer  (inviscid)  solutions  by 
Heisenberg  (Ref.  52) .  The  classical  ^  is  uniformly  slowly  varying, 
whereas  the  ^  given  here  is  moderately  varying  fory'~  0(1),  and  is  a 
rapidly  varying  function  of  the  physical  coordinate  y  This 
solution  tends  only  asyirptotically  for  9  »  1  to  a  slowly  varying 
function.  ' 


Also  we  note  that  solution  (A, 31)  coincides  with  the  leading  term 
of  the  classical  outer  solution,  which  further  substantiates  the 
asymptotic  method  applied  here. 


Solution  <1^ 

Repeated  use  of  the  theorem  assures  that  the  third  solution 
^  satisfies  a  first  order  linear  differential  equation  to  be  deduced 
from  (A. 13) .  This  assertion  is  realized  by  the  method  of  variation  of 
constants 


(A.33) 


where  ^ and  ^^are  functions  of  f  to  be  determined.  Having  eliminated 
4)^  between  the  original  equation  (A. 13)  and  the  supplementary  condition 


'5‘S 


which  is  standard  to  this  method,  we  are  lead  to  the  equation  for 


n.  ■  >f  4  -  4 


This  equation  is  integrated  twice 
gives  ^ 

<5  = 


,  and  its  substitution 

a 


into 


(A.33) 

(A. 34) 
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where  W  is  the  V/ronskian  formed  by  and 


yr- 


(A. 35) 


It  is  easily  confirmed  from  solution  (A. 34)  for  that  it  has  an 
asymptotic  form  of  exponential  growth  with  ^  ,  In  fact,  using  the 
same  approximation  as  has  been  used  previously  in  deriving  (A. 31) ,  we 
have  for  (A. 34] 


with 


/ 


-sU 


ZS^^,(U-c) 


(A. 36] 


In  deriving  the  second  row,  the  asymptotic  expressions  (A. 24)  and 
(A. 31)  for  and  <8^  have  been  used.  Relationship  (A. 36)  is  again  in 
qualitative  agreement  with  its  classical  equivalent,  as  is  the  case 
for  the  decaying  solution  .  The  classical  equivalent  is,  from 
Ref,  51,  „  ^ 


CJ  -CO  -et 


■/3 


(A. 37) 


with  the  same  nomenclature  as  before  with  the  Hankel  function  of  the 
second  kind  ,  The  asymptotic  expression  of  (A.25)  for  7'»  1  is 


(A.38) 


in  agreement  with  (A. 36)  with  limited  accuracy  in  the  velocity 
profile.  ^ 


The  set  of  solutions  which  has  been  constructed  has  a 
noteworthy  characteristic  for  the  Wronskian  that  simplifies  the 
analyses  which  follow 


This  formula 


K 

% 

is  easily  checked  by  noting  the  following 


(A. 39) 

relationship 

>y 
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which  holds  for  J-  0,  If  2. 


Solution 


Since  all  homogeneous  solutions  are  exhausted,  the  fourth 
solution  of  equation  (A. 12)  is  the  particular  solution  of  the 
equation  with  C,  =  1 


-  c( U- cj  ^  ^ 


(A.40) 


The  method  of  variation  of  constants  is  used  here  also,  since  all  of 
the  homogeneous  solutions  through  are  available.  The  theorem 
states  that  no  differential  equation  needs  to  solved.  A 
manipulation  using  the  key  property  (A. 39)  leads  to  the  solution 


^  r  \ 

with  -5  =  V 

where  9  =7»  1  is  a  point  beyond  which  the  asyn^stotic  expressions  for 
^  through are  valid.  Taking  the  derivatives  successively,  we  have 


(A. 41) 


(A.42) 


r  \  ^  r  ' 


r 


7 


(A.43) 


This  result  can  be  checked  by  using  (A, 38) ,  (A, 32) ,  and  (A. 33) , 

The  function  in  expression  (A. 41)  behaves  as  a  slowly  varying 
function  for  f  »  \  under  the  same  conditions  as  invoked  in  deriving 
the  asyirptotic  expression  {A. 31)  for  <P^  •  The  resulting  expression  is 


f  ^2 

<^^~-(u-c)'\{D-^)  (^>7»i) 

For  comparison,  the  classical  equivalent  to  this  function  is 

(o-c) 


(A.44) 


y  -z 

^  J-  ^  (A. 45) 

where/ is  defined  in  (A.24) .  The  leading  terms  of  the  asymptotic 
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expression  agree  with  the  leading  terms  of  the  classical  outer 
solution. 


Because  of  (A. 43) t  we  can  show  the  following  relationship  for  the 
Wtonskian  formed  by  the  four  solutions  ^through 


= 

■  * 

■  w 

r«* 

; 

V  *  • 

(A. 46) 


Order  Approximation 


The  full  solution  of  the  Orr-Sommerfeld  equation  is  constructed 
from  the  series  (A. 6)  by  successive  approximation  starting  with 
equation  {A.13)  with  solution <^’known.  In  practice,  the  nth-order 
correction  is  obtained  by  solving  the  inhomogeneous  equations 
(A.IO)  following  the  same  procedure  as  has  lead  from  equation  (A. 40) 
to  solution  (A. 41).  The  result  is 


with  ~ 


<'*'0  .  \ 


-\ 


{A. 47) 
(A.48) 


9  ^  ^ 

Using  these  formula,  the  same  approximation  as  used  in  deriving 
the  asynptotic  expressions  for  ^  and  ^  ,  we  get 

(for^>/»l)  (A.49) 

This  expression,  when  substituted  into  A. 5,  with  ^(/;1,2)  given 
by  (A, 44)  and  (A. 31)  respectively,  coincides  with  Heisenberg's 
inviscid  outer  solutions  (Ref.  52).  The  series^  solutions  for  F,  and 
liare  uniformly  convergent^if  the  lower  bound  y=e/is  taken  such  that 
U-c  varies  only  slowly  for  y  <  y.  In  fact,  then,  the  series  sum  up  to 
yield  . 


=  (c^-c)  roj>A  ^ 


(A.50) 
(A. 51) 


In  closing,  this  solution  of  the 
Tollmien  wave  has  been  calculated 
(Ref.  2) . 


Orr-Sommerfeld  equation  for  the 
and  verified  by  Tsuge  and  Sakai 
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^pendix  B 

A  sununary  of  basic  characteristics  of  five  forms  of  2-D  spatial 
oscillations  in  a  parallel-flow,  Blasius  boundary  layer,  viscous,  flat 
plate,  with  0<y<ao 


1.  The  Tollmien-Schlichtinj  instability  Efiye 
is  the  fundamental  wave  in  a  finite  set  of 
elgennodes.  The  T-S  wave  is  the  only  known 
unstable  eigenwave  for  an  incompressible,  flat 
plate  boundary  layer.  The  higher  modes  are 
all  heavily  damped  for  the  Blasius  layer.  The 
main  physical  features  Include  two  viscous 
layers  (a  viscous  sublayer  near  the  wall,  and 
a  critical  layer  near  where  U(y^  )=^Real(the 
phase  speed)).  The  T-S  wave  propagates 
downstream  and  may  grow  or  decay  depending  on 
the  boundary  layer,  Reynolds  number,  and 
frequency.  Across  the  sublayer,  the  Reynolds 
stress  jumps  from  zero  to  some  value.  Across 
the  critical  layer,  that  stress  jumps  back  to 
zero.  This  Reynolds  stress,  if  viscous 
dampening  is  not  too  large,  can  lead  to  growth 
of  the  wave. 


Figure  B,1 


2.  The  decavinz  standing  wave  decays 
exponentially  in  the  streamwise  direction  and 
oscillates  sinusoidally  in  time  as 
v=f(y)exp(-px-itft) .  This  oscillation  does  not 
propagate  in  the  streamwise  direction!  the 
phase  speed  is  pure  imaginary.  Above  the 
boundary  layer,  the  standing  wave  behaves  as 

v=[Aexp(-my)+Bsinpy+cospy]exp(-f,x-ia/t) . 

Far-away  from  the  boundary  layer,  this 
fluctuation  is  irrotational .  The  limit  of  a 
standing  wave  as  is  a  Stokes  wave 
u(y,t)=f(y)exp(-l<At) ,  v=0,  for  an  oscillating 
freestream.  The  “collision"  of  vortical  and 
other  freestream  disturbances  with  the  leading 
edge  can  excite  a  spectrum  of  these  decaying 
waves.  Interactions  between  the  trailing  edge 
and  disturbances  are  also  believed  to  excite 
these  waves.  Steady  waves  of  this  form  are 
also  Induced  by  flow  over  a  wavy  wall  with  a 
leading  edge. 


Figure  B.2 
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3.  The  growing  ttendlng  wave  it  the 
exponentially  growing  counterpart  to  (2).  If 
invitcid,  oacillationa  (2)  and  (3)  are  related 
by  Although  the  viscoua  terms  do  not 
transform  this  way,  calculations  show  that 
this  relationship  is  a  good  guide  for  large 
Reynolds  numbers.  The  steady  limit  of  this 
wave  behaves  fundamentally  different  than  (2). 
The  phase  speed  it  pure  imaginary.  This 
oscillation  it  believed  to  be  excited  by 
leading  and  trailing  edges.  Like  (S).  this 
oscillation  represents  an  upstream  influence 
of  downstream  b.c.  in  numerical  calculations. 
Wall  wavlnest  can  also  excite  these 
fluctuations . 


Figure  B.3 


4.  Decaying  vortical  fluctuations  propagate 
downstream  with  a  speed  slightly  greater  than 
the  freestream.  and  decay  slowly  as 
v*^(y)oxp(i^x-c^t))exp(^x)  where 

+4fe*) /Rj+. . .  A  passive  viscous  sublayer 
forms  at  the  wall.  The  smallness  of  the 
velocity  fluctuations  in  the  boundary  layer  is 
caused  by  the  formation  of  a  layer  of 
vorticity  at  the  “edge*'  of  the  boundary  layer 
which  induces  a  flow  in  the  direction  opposite 
to  that  induced  by  the  vortical  freestream 
disturbances.  Outside  of  the  boundary  layer, 
the  flow  behaves  as 

v'°[Aexp(-my)-EBsin^y-*-cospy]exp[i><(x-ct)] .  Far 
above  the  boundary  layer,  the  fluctuations  are 
rotational .  in  contrast  to  waves  (2)  and  (3) . 


Figure  B.4 


5.  The  traveling  vortical  wave  grows 
explosively  like  expC-t-R^x).  It  propagates 
upstream  at  a  speed  approximately  equal  to  -D^, 


and  represents  upstream  diffusion  of 
vorticity.  Difficult  to  calculate  due  to  the 
very  high  frequency  oscillation  exp(iR^y)  in 
the  boundary  layer.  This  oscillation,  like 
waves  _(2)-(4),  also  exists  for  a  uniform  mean 
flow,  0~1 ,  in  contrast  to  the  stability  waves 
for  which  a  boundary  layer  must  be  present. 
Above  the  boundary  layer,  this  oscillation 
behaves  as 
v=[Aexp(-iBy)+Bcos^y+sin^]exp[io<(x-ct)]  .  This 
wave  is  one  of  the  upstream  influences  of  a 
downstream  boundary  condition  in  a 
calculations!  domain.  It  is  one  of  the 
upstream  Influences  of  leading  and  trailing 
edges,  and  is  perhaps  an  additional  viscous 
diffusion  effect  in  non-parallel  flows. 


Figure  B.5 
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NOMENCLATURE 


English 
a,  b 


C=C,.  +1C{ 

D=d/dy 

f(y)=f^+ifi 

F  (a)/  F  (a) 

h> 

l'”,  l'" 
p 

R»  H 

s=i'»< 
s 


i3'  H 


& 

P 

t 

U,  V,  w 

U{y) 

V,  w 


w 

X 

y 


y 

z 


convergence  parameters  in  generalized  Fourier 
transforms,  equation  (2.6) 
defined  in  equation  (A. 24) 

constants  in  equation  (3,19a) 
complex  phase  speed 
complex  constants  in  equation  (3.B) 
ordinary  derivative 

complex  amplitude  of  the  longitudinal  velocity 
forcing  function 

Fourier  coefficients  of  the  forcing  function 
defined  in  equation  (4.20) 
defined  in  equation  (A. 2 8) 

constants  in  the  partial  fractions  of  equation  (3.21) 
(-1)''“' 

contribution  from  the  jth  pole  as  defined  in 
equation  (4.2) 

operators  defined  in  eqns.  (A. 11, A. 12) 
disturbance  pressure 

Reynolds  number  based  on  characteristic  thickness 
of  the  boundary  layer 
Laplace  parameter 

real  value  of  s  lying  to  the  right  of  all  poles 
values  of  s  defined  in  equation  (4,11) 
transformed  variable  in  equation  (A, 28) 
time 

disturbance  velocities  in  the  x,  y  and  z  directions 
mean  velocity  in  the  streamwise  direction 
mean  x-velocity  in  the  freestream 
roeein  velocities  in  the  y  and  z  directions; 
see  equation  (A. 17)  for  use  of  V,  W  as 
transformed  variables  in  Appendix  A 
Hronskian 

coordinate  parallel  to  plate  in  the  streamwise  direction 

coordinate  normal  to  the  plate 

y-value  of  the  boundary  layer  edge 

y-value  for  the  asymptotic  matching 

spanwise  coordinate 


Greek  Script 
of  X  -wavenumbe  r 

X -wavenumbers  of  the  discrete  modes 
^  y- wavenumber 

r  z-wavenumber ;  contour  of  integration 

V*  Laplacian  operator 

S  characteristic  thickness  of  the  boundary  layer; 

delta  function 
cofactor 

A  defined  in  equation  (3.18) 

f=l/R  inverse  of  Reynolds  number 
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Z  j/2 

[«^  +iR(o<-afl  ]  •  ;  see  equation  (A.lSarb) 

for  definition  in  Appendix  A 
kinematic  viscosity 

disturbance  vorticity  in  the  z-direction 
complex  amplitude  of  the  normal  velocity  disturbance 
four  independent  solutions  of  the  Orr-Sommerfeld  equation 
defined  in  equation  (3.17) 
integration  contour 

-•»<>»<«»,  imaginary  part  of  s  lying  to  the  right 
of  all  poles;  variable  defined  in  equation  (A, 13) 
frequency 

vector  defined  in  equation  (4.13) 
defined  in  equation  (4.18) 

Superscripts.  Subscripts.  dM  Miscellaneous  Notation 
(  )  time  average 

e  boundary  layer  edge 

p  particular  integral 

r,  i  real  and  imaginary  parts  of  a  complex  variable 

()^  partial  derivative  with  respect  to  x 

^  Fourier  transform 

(0)  evaluated  at  x=0 

c.l.  classical  solution 


A 


y)-^+i;^- 

,i=l,2,3,4 

(0) 
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